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We present a theoretical study ol recent laser-alignment and mixed-field-orientation experiments 
of asymmetric top molecules. In these experiments, pendular states were created using linearly 
polarized strong ac electric fields from pulsed lasers in combination with wealt electrostatic fields. 
We compare the outcome of our calculations with experimental results obtained for the prototypical 
large molecule benzonitrile (C7H5N) [J. L. Hansen et al, Fhys. Rev. A, 83, 023406 (2011)] and 
explore the directional properties of the molecular ensemble for several field configurations, i. e., for 
various field strengths and angles between ac and dc fields. For perpendicular fields one obtains 
pure alignment, which is well reproduced by the simulations. For tilted fields, we show that a fully 
adiabatic description of the process does not reproduce the experimentally observed orientation, 
and it is mandatory to use a diabatic model for population transfer between rotational states. We 
develop such a model and compare its outcome to the experimental data confirming the importance 
of non-adiabatic processes in the field-dressed molecular dynamics. 



I. INTRODUCTION 

Controlling molecular motions has direct impact in 
a wide variety of molecular sciences, including stereo- 
chemistry molecular- frame investigations of geo- 
metric and electronic properties, such as photoclcctron 
angular distributions [5H7j and high-harmonic genera- 
tion [1, Q, as well as for diffractivc imaging of gas- 
phase molecules [l^, [HI , aiming at recording the "molec- 
ular movie" [l^l- Recently, there has been tremendous 
progress in the control of the translational [l^ and ro- 
tational motions of even complex molecules. For 
extremely well controlled ultracold alkali dimers, direct 
quantum effects in the stcrcodynamics of molecular reac- 
tions have recently been observed (Tsj . 

Angular confinement of molecular ensembles is referred 
to as alignment ~ the confinement of molecule-fixed axes 
along laboratory-fixed axes - and orientation - adding 
a well-defined direction. Traditionally, these two lev- 
els of angular control have been separated: strong ac 
fields from pulsed lasers have been used to create align- 
ment [l3|, whereas state-selection [l^ and brute- force 
orientation using strong dc electric fields [l^, [2l| have 
been used to create orientation (which typically also cre- 
ates alignment). About a decade ago it was realized that 
strong simultaneous alignment and orientation could be 
created using combined ac and dc electric fields 2^ 23 . 
This has been experimentally verified in a few cases 24- 



|26[. Adding methods to control the translational mo- 
tion one can perform a quantum-state selection before 
the alignment and orientation experiment [l^ [13, [28| . 
This two-step approach has allowed the creation of un- 
precedented degrees of one-dimensional (ID) and three- 
dimensional (3D) alignment and orientation even for 
complex asymmetric top molecules [l^, [29I - I3H . 

Theoretical studies of the rotational spectra in the 
presence of combined electrostatic and non-resonant ra- 
diative fields have been restricted to linear and symmet- 
ric top molecules so far [2^ [s^ • Recently, some of 
us developed the theory for asymmetric top molecules 
in combined fields [s^l- When the static and linearly 
polarized laser field directions arc tilted with respect to 
one another, the symmetries of the corresponding Hamil- 
tonian are significantly reduced. This is the most gen- 
eral field configuration, its theoretical treatment being 
most challenging. Since each external field interaction 
breaks different symmetries of the field-free Hamiltonian, 
the order in which the fields are turned on determines 
the evolution of the field-dressed rotational dynamics. 
The labels of the eigenstates at a certain field configura- 
tion obtained by adiabatic following depend on the path 
through parameter space, i.e., monodromy is observed. 
For linear and symmetric top molecules exposed to ex- 
ternal fields, the phenomena of classical and quantum 
monodromy have been encountered in the corresponding 
dressed spectra (s^ - fSGj . The labeling procedure is also 
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numerically very demanding due to the large amount of 
genuine and avoided crossings occurring between adja- 
cent states as one of the field parameters is varied. 

In the present work, we describe the recent experimen- 
tal results obtained for the alignment and orientation of 
asy mmetric top molecules in combined fields 
l3l| . In the framework of the rigid rotor approxima- 
tion, we perform a full investigation of the field-dressed 
eigenstates for those field configurations considered in 
the experiment. The Schrodinger equation for the ro- 
tational/pendular states is solved separately for each ir- 
reducible representation by expanding the wave function 
in a basis with the correct symmetries. Our theoret- 
ical analysis also includes i) a diabatic approximation 
to account for population transfer through the avoided 
crossings as the laser intensity is varied; ii) the velocity 
distribution of the ions after the Coulomb explosion; and 
iii) a volume effect model to describe the fact that not 
all the molecules feel the same laser intensity because the 
orienting and the detection laser pulses have finite spa- 
tial intensity profiles. The field-dressed eigenfunctions 
are weighted with the known relative state populations 
in the molecular beam (29| . Then, we compute the an- 
gular probability density functions for different field con- 
figurations and their velocity- mapping images (VMIs). 
The recently performed alignment and orientation exper- 
iment for benzonitrile (BN, C7II5N) molecules [l,[3l| pro- 
vides us with experimental data that is very well suited 
to present and discuss our theoretical model in a com- 
parative study. We numerically compute the alignment 
and orientation for an ensemble of quantum-state se- 
lected benzonitrile molecules, and compare our results 
to the experimental data. For comparison with a cold 
thermal ensemble - without state selection ~ we also pro- 
vide the corresponding results for a benzonitrile sample 
at 1 K. In particular, for perpendicular fields, we obtain 
good agreement between the computational and exper- 
imental results for the degree of alignment. For tilted 
fields, we show that a fully adiabatic description of the 
rotational/pendular dynamics cannot reproduce the ex- 
perimental results for the mixed-field orientation of ben- 
zonitrile. When a diabatic model is implemented for the 
treatment of the avoided crossings, our theoretical study 
reproduces with reasonable accuracy the experimental 
degree of orientation. Hence, we demonstrate the impact 
of non-adiabatic processes on the field-dressed molecu- 
lar dynamics. We have developed a general theoretical 
description of alignment and mixed-field orientation for 
asymmetric tops in long pulses of strong ac electric and 
weak dc fields. 

The paper is organized as follows: In Sec. |ll]thc rel- 
evant experimental details are described. The theoret- 
ical model is presented in Sec. IIIIl which includes the 
discussion of the rigid rotor Hamiltonian and its symme- 
tries, the diabatic model to treat the avoided crossings, 
the screen projection of the 3D probability densities and 
the experimental observables. The theoretical results for 
the alignment and orientation of a beam of benzonitrile 



molecules are compared to the experimental data in Sec. 
IIVI The conclusions and outlook are provided in Sec. |Vl 

II. EXPERIMENTAL DETAILS 

A detailed description of the alignment and mixed- 
field orientation experiments is given elsewhere [1, [2^, 
Isij . Briefly, a pulsed, cold molecular beam of benzoni- 
trile molecules seeded in helium is expanded from an 
Even-Lavie valve into vacuum. The molecular beam is 
skimmed before entering a 15 cm long electrostatic de- 
flector 41 cm downstream from the nozzle. The deflector 
disperses the molecules in the beam according to their ef- 
fective dipole moments, creating a vertically varying dis- 
tribution of quantum states in the probe region 77 cm 
downstream from the nozzle. Alignment and orienta- 
tion is induced by the dc electric held of the VMI spec- 
trometer and by a strong Nd:YAG laser pulse (10 ns, 
1064 nni) and probed using ion imaging of CN+ frag- 
ments following Coulomb explosion of benzonitrile with 
a strong ultrashort Ti:Sapphire laser pulse (30 fs, 800 nm, 
5.4 X 10^^ W/cm^). Alignment and orientation experi- 
ments were performed for the undeflected beam and for 
the quantum-state selected sample at a vertical height of 
1.75 mm; see Figure 2 of reference [Sll for details. 

III. THEORETICAL MODEL 

An exact theoretical quantum description of the align- 
ment and orientation experiments is very demanding, 
since it requires the solution of the corresponding time- 
dependent Schrodinger equation for each state populated 
in the molecular ensemble. Instead, we retreat to a quasi- 
static description parametric in the fleld strength and an- 
gles, and solve the time- independent Schrodinger equa- 
tion for several field configurations. A diabatic model, 
based on the field-free symmetries, is used to account 
for the population transfer through the avoided crossings 
encountered as the YAG pulse intensity is varied. 

A. The Hamiltonian 

We consider a nonresonant laser field of intensity / 
and linearly polarized along the Z^-axis of the laboratory 
fixed frame (LFF), {Xl,Yl, Zl), and an homogeneous 
electrostatic field of strength Es contained in the XlZl- 
plane and forming an angle f3 with Z^. The rigid rotor 
Hamiltonian of a polar asymmetric top molecule exposed 
to this field configuration is given by 

H = JxmBxm + J'LBy,, + Jz^Bz,, - Esfi cos 03 
2ttI 

{a^"'^'' cos^ e + a^""^^' sin^ 9 sin^ x), (1) 

c 

with Bxm, Bymi and Bzm being the rotational con- 
stants. The molecule or body fixed frame (MFF) 
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{XmjYm, Zm) is defined so that the permanent electric 
dipole moment /x is parallel to the Zjv/-axis. The LFF and 
the MFF are related by the Eulcr angles (0, 8, y)^i\. 
The polarizability tensor is diagonal in the MFF with 
components an with i = Xm, Ym, Zm, and the inter- 
action with the laser field depends on the polarizability 



the Hamiltonian (1) is invariant under C^", a rotation of 



anisotropics 



OLi. 



XmiYm, Z 



M- 



The 



angle between the static electric field and the molecular 
^M-axis is 9s with cos 9s = cos (3 cos ^ + sin ^ sin 9 cos 0. 
This study is restricted to molecules with the permanent 
dipole moment parallel to one of the principal axis of 
inertia, as it is the case for benzonitrile. 

Let us shortly summarize the approximations and as- 



sumptions made in the derivation of the Hamiltonian (1~ 
We perform a non-relativistic description within the 
framework of the rigid rotor approximation, assuming 
that the electronic and vibrational structures are not af- 
fected by the external fields. In addition, we presuppose 
that the laser is non-resonant and that the inverse of the 
oscillation frequency is much larger than the rotational 
period of the molecular system and pulse duration. Thus, 
we can average over the rapid oscillations, so that the in- 
teraction of this field with the molecular dipole moment 
is zero, and only the interaction with the polarizability is 
left. Additionally, since the YAG laser pulse duration is 
much larger than the time scale of the rotational dynam- 
ics, we assume that the alignment and orientation pro- 
cesses are, in principle, adiabatic, and we take the time 
profile of the pulse as a constant equal to 1 . The validity 
of this assumption will be discussed in detail through the 
paper. The spatial dependence of the laser intensity will 
be also taken into account by means of a volume effect 
model, see below. 

The configuration of the fields determines the symme- 
tries of the Hamiltonian 1(1) For the field- free case, they 



are given by the spatial group S0(3) and the molecu- 
lar point group D2 - i.e., the Fourgroup V - consisting 
of the identity and the two-fold rotations C\ around the 
MFF i-axis, with i = Xm, Ym and Zm- The Schrodinger 
equation associated with the Hamiltonian (1) can not 



be solved analytically, and only the total angular mo- 
mentum, J, and its projection onto the Z^-axis of the 
laboratory fixed frame, M , are good quantum number, 
whereas if, the projection of J on the ^M-axis of the 
molecular fixed frame, is not conserved. Since an external 
field defines a preferred direction in space, the symme- 
tries of the corresponding Hamiltonian are reduced com- 
pared to the field-free case. Here we consider three field- 
configurations. For a single static field parallel to the 
LFF- Z^-axis, the symmetry operations arc , any 
arbitrary rotation around the Z^-axis and a reflection 
in any plane containing Z^, thus, M is still conserved. 
For a certain |M|, there are 4 irreducible representation, 
and for Af 7^ the states with M and —M are degener- 
ate. If the molecule is exposed to both fields, with the 
electric field being rotated away from the Z^-axis, the 
azimuthal symmetry is lost and M ceases to be a good 
quantum number. In the perpendicular case, /3 = 90°, 



TT around the X^-axis, and the refiection aX[^ZL {X^Z^ 
is the plane containing the fields); in consequence there 
are 8 irreducible representations. For tilted fields with 
(3 7^ 90°, two symmetries are left, txlZl and C^"^, and 
the corresponding group has only 4 irreducible represen- 
tations. For a detailed description and analysis of the 
symmetries for all possible situations, we refer the reader 
to the recent works [ssi |38|. 

The field-dressed eigenstates of an asymmetric top are 
characterized by the avoided crossings appearing between 
levels of the same symmetry as one of the parameters of 
the field configuration, i.e., Es, I, or /?, is varied. For 
non-coUinear fields, the small amount of irreducible rep- 
resentations implies an eigenstate diagram with a high 
degree of complexity due to the large number of avoided 
crossings. Furthermore, a small degree of asymmetry on 
the inertia tensor facilitates the appearance of avoided 
crossings on the corresponding field-dressed states [39j . 
For the correct analysis of the computational results, 
these avoided crossings should be distinguished from gen- 
uine ones taking place between levels of different symme- 
tries. Hence, the time- independent Schrodinger equation 
associated to the Hamiltonian (1) is solved for each irre- 
ducible representation by expanding the wave function in 
a basis that respects the corresponding symmetries [s^. 

The field-free states are identified by the quantum la- 
bels Jk^k^M, with J and M being good quantum num- 
bers, and Ka and Kc the projections of J onto the Zm- 
axis of the molecular fixed frame in the oblate and prolate 
limiting cases [i^l, respectively. For reasons of address- 
ability, we will denote the field dressed states by means 
of these field-free labels indicating if an adiabatic or dia- 
batic picture has been used. 

To illustrate the molecular dynamics in the state- 
selection process using the electric deflector, we plot 
the Stark energies and expectation value {cos 9) = 
{JKaK^M\ cos9\JKa^K^M) of the populated states of ben- 
zonitrile in |FigurcT] (a) and (b) as a function of the elec- 
trostatic field strength. This plot covers the range of 
electric field strengths present in the deflector. It includes 
84 individual rotational states accounting for 98 % of the 
population of the molecular ensemble. This level dia- 
gram shows a complex structure with both genuine and 
avoided crossings. In the weak electrostatic fleld regime, 
both high- and low-fleld seekers are encountered, whereas 
in the pcndular limit, and, in particular, in the regime of 
interest for the state selection experiment, all these lev- 
els are high-field-seekers. For Es > 50 kV cm~^, most 
of these populated states present a significant orientation 
with {cos9) > 0.5, and 32% of them have a strong one 
with (cos 61) > 0.8. The avoided crossings have a strong 
impact on the character of the involved states, and (cos 9) 
suffers large variations over tiny ranges of electrostatic 
field strength, sec |Figure l] (b). Note that the dc field 
applied in the orientation experiment, 286 Vcm^^. is 
very small on the scale of this figure. 
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Fig. 1: (a) Energies and (b) expectation value {cos 9) as a 
function of the static field strength of the populated rotational 
states for benzonitrile. 



B. The diabatic model 

After entering the region of the extractor field, i.e., the 
static electric field of the VMI spetrometer, the molecules 
are exposed to a laser field with increasing intensity 
(YAG laser pulse). As the laser intensity varies, a certain 
state may undergo several avoided crossings with levels of 
the same symmetry. The presence of these avoided cross- 
ings as well as their diabatic or adiabatic nature have a 
strong impact on the outcome of the experiment. 

Several theoretical studies have analyzed in detail the 
character of such avoided crossings for molecules exposed 
to an electrostatic field by using different adiabaticity 
criteria (sol . l4l| - |43| . Their main and common conclusion 
is that the assumption of a fully adiabatic dynamics is 
incorrect as no general statement can be made about 
the character of the avoided crossings. These works sug- 
gest that an investigation of the character of the avoided 
crossings encountered as / is varied should be mandatory 
for a correct description of the experimental results. To 
do so, we use the following adiabatic passage criterion [iJl 





dH' 




{• 


at 





< 1, 



(2) 



where E^ and Ej arc the eigenenergies for the states i and 
J and H' is the interaction term. Due to the large amount 
of avoided crossings in the laser field-dressed states, a 



systematic and detailed study of the adiabatic or diabatic 
character of all of them is unfeasible. Thus, we employ a 
simple diabatic model that provides an approximation to 
the dynamics, and determines the population transfer as 
the field parameters are varied. Note, that some previous 
works have proposed different population transfer models 
based on symmetry considerations [4l|, IH, |45| . 

If the fields are not collinear and the polarization of 
the laser pulse is parallel to the Z^-axis, the electrostatic 
field induces the coupling of states with different field-free 
values of the quantum number A/. For the experimen- 
tally used electrostatic field strength, Es ~ 286 V cm~^, 
the Stark interaction is much weaker than the laser field 
one. Thus, the hybridization of the quantum number 
M is so small that for a certain level (M^) remains al- 
most unperturbed and equal to its field free value. Our 
diabatic model consists in assuming that i) an avoided 
crossing between two levels with different field-free val- 
ues of M is considered as being crossed diabatically; and 
ii) crossings between levels with the same field-free value 
of M arc passed adiabatically. To illustrate the valid- 
ity of this diabatic model we have analyzed two exem- 
plary avoided crossings of the field-dressed spectrum with 
/? = 45° and Es = 286 V cm"-'^ by means of the adiabatic 
passage parameter 77 of equation (2)| In this case, the 



interaction term is given by H' ~ (27r//c)(a 



cos^ I 



6 sin x). For the time profile of the YAG pulse. 



a sm 

we use /(i) = /exp(-tV(2cr2)) with / ^ 5x10" Wcm"^ 
and a = 4.25 ns (FWHM=10 ns). [Figure 2| (a), (d), 
(b), (e) and (c), (f) depict the energy, (cos^ 6*) , and ry 
as a function of / for the avoided crossings between the 
states 202O and 2o2l, and 4o44 and 821 1, respectively. 
The levels 202O and 2o2l suffer an avoided crossing for 
/ « 3.5 X 10^ W cm-^ with /\E = 5.5 x 10"^ cm^i 
and rj = 192. Hence, we conclude that it is crossed dia- 
batically. In contrast, for the avoided crossing among the 
states 4o44 and 32il occurring at / w 3.82x 10^° W cm"^, 
we obtain AE" = 8.69 x 10"^ cm"^ and 77 = 0.12, which 



according to equation (2)| is an intermediate case - nei- 
ther diabatic nor adiabatic. Using a fully adiabatic pic- 
ture these two states are labeled 4o44 and 821 1, but using 
the diabatic model their labels are different since (M^) = 
3.999 and 4.007, respectively, for / = 3.6 x 10^° W cm-^, 
which explains the mixing between both levels. The cri- 
terion 1(2)1 or the Landau- Zencr formula do not classify 
this avoided crossing as being diabatic or adiabatic, but 
within our approximation we consider it to be adiabatic. 

For tilted fields (/3 7^ 90°), a certain avoided cross- 
ings may involve two states one being oriented and the 
other one antiorientcd. Hence, considering it as diabatic 
or adiabatic has significant consequences on the final re- 
sult for orientation. Of course, our diabatic model is an 
approximation, since many avoided crossings are encoun- 
tered that do not fall clearly into the class of adiabatic 
or diabatic crossings according to our adiabatic passage 
criterion. Generally, these cases would require to solve 
the time-dependent Schrodinger equation including the 
time profile of the YAG pulse. 
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Fig. 2: Energy E ((a) and (d)), expectation value (cos '^ 6) 
((b) and (e)), adiabatic passage criterion parameter (2) ((c) 
and (f)) for the levels 2o2 (dash-dotted line) and 2o2l (dashed 
line) states; and for the states 4o44 (dash-dotted line) and 821 1 
(dashed line) states as a function of the YAG laser intensity, 
for Es = 286 V cm~^ and ^ = 45°. See the text for more 
details. 



If the fields are perpendicular, the corresponding 
Hamiltonian (1) has 8 irreducible representations, and 



the amount of avoided crossings is significantly reduced 
compared to the /3 ^ 90° configuration. For the pop- 
ulated states of one irreducible respresentation, the ex- 
pectation value (cos^ 9) is plotted on [Figure 3| as a func- 
tion of the laser intensity, for Es = 286 V cm~^ and 
/3 = 90°. The inset of this figure shows the energy levels 
in the strong-laser-field regime. For the laser intensities 
at which the experiment is performed, most of the popu- 
lated levels are characterized by a pronounced alignment. 
In particular, we find (cos^ 9) > 0.8 for all these states 
and / > 4.1 X 10^^ W cm^^. Since, the character - re- 
garding alignment- of two states that had previously suf- 
fered an avoided crossing is very similar, passing through 
it diabatically or adiabatically is not crucial for the final 
alignment result. 
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Fig. 3: Expectation value (cos^ 9) and energies (inset panel) 
as a function of the YAG pulse intensity of the populated 
rotational states for benzonitrile, for Es = 286 V cm^^ and 
/3 = 90°. Only one irreducible representation is presented. 



The projection of the probability density on the 
screen detector 



When a molecule is multiply ionized using an intense 
ultrashort laser pulse, it dissociates due to Coulomb re- 
pulsion and the created ionic fragments are collected in a 
2D space resolving detector. Within the axial recoil ap- 
proximation the momentum of the CN+ fragments cre- 
ated in this process is parallel to the C-CN bond di- 
rection before ionization. In order to calculate the screen 
image from the molecular wave functions we define a new 
reference frame {x,y,z), containing the detector screen 
on its 2/z-planc. This frame is obtained by rotating the 
LFF by an angle 90° — /3 about the l^-axis. Note that 
the electrostatic field is parallel to the x-axis, and the 
polarization vector of the linearly polarized YAG laser 
lies in the a;z-plane, forming an angle /3 with the x-axis. 

An asymmetric top molecule in a certain state, char- 
acterized by the wavefunction 'i-y(9,(j),x)i with 7 = 
JxaK^M, is traveling towards a screen (velocity in the 
x-direction). Since the direction of the molecular dipole 
moment is independent of the angle x, the probability 
for the molecular ^M-axis to be oriented according to 
(6*, 0) is given by integrating in x the absolute square of 
the wave function. The angular distribution p^{9, (f))dQ 
provides a measure of the amount of ions ejected into 
the solid angle dil = sin 9 d9 dcf>. It is related to a spa- 
tial distribution on the 2D screen by p^{y^z)dydz = 
p~f{9{y,z),(j){y,z))\J\s\n-{9{y,z))dydz, with J being the 
Jacobian of the transformation between the coordinates 
(6', 0) and {y,z). Assuming that all the ions have the 
same velocity, v in absolute value, this transformation 
reads as 



y ^ a sm 9 sm cp 

z — a(cos 9 sin /3 + sin 9 cos 4> cos / 



(3) 



where a = vtf, and tf is the time of fiight to reach 
the screen. Note that with the definition |(3)| for the y- 
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coordinate, any point (y, z) corresponds to two different 
orientations {9, (p) and {6, 180° — </>). Thus, the screen im- 
age is the sum of two projection images, each obtained 
for a restricted range of 0. 

An important ingredient, that should be taken into ac- 
count to obtain reahstic screen images, is the ahgnment 
selectivity of the probe laser. We approximate the ef- 
fectivity of the dissociation process by the factor cos^ a, 
with a being the angle between probe laser polarization 
and the C-CN bond direction [301 ■ Note that this ap- 
proximation for modeling the probe selectivity as cos^ a 
is in accord with the experimental observations. Thus, 
the 2D screen spatial distribution is given by 



Pi{G{y.z),(t){y,z)) 



(4) 



where Ai(y, z) is the mentioned aligment selectivity fac- 
tor, and the index i indicates the probe pulse polar- 
ization. For a linearly polarized probe beam in the 
x-direction, we obtain Ai(y^z) = 1 — y^/a^ — ^^/i^, 
and in the z-direction Ap(y, z) = z/a. For tilted fields 
(/3 7^ 90°), a circularly polarized probe in the xz-plane, 
ensures that any molecule is ionized and detected with 
the same probability independently of j3, and it gives 
Ac(?/, z) = 1 — y^ja}. The apparent singularity of 
equation (4) at + z^ = a} will disappear when we 



integrate over different ion recoil velocities. 

So far we have assumed, that all ions acquire the same 
recoil velocity. Experimentally, however, they follow a 
certain distribution D{a), which is assumed to have only 
nonzero values at positive velocities. The screen image is 
obtained by averaging over all these velocities with their 
corresponding weights as 



pUy, z; a)D{a) da. 



with i = l,p and c depending on the probe polarization. 
The experimental velocity distribution D{a) for the CN+ 
ions created in the Coulomb explosion of benzonitrile is 
shown in [Figure 4"! The second and third peak are the 
two Coulomb explosion channels relevant for determin- 
ing the orientation of the C-CN axis. The two peaks are 
fitted to a combination of two Gaussian functions, after 
subtraction of background, and used further on in the 
model. Since a is not the velocity but a distance propor- 
tional to it, we have rescaled the abscissa such that the 
maximum of D{a) is at a = 1. 

In the experiment, the molecular beam contains 
molecules in different quantum states. A certain level 
7 = Jk^k^M has a relative weight = Wj^^^^^^m 
within this ensemble, which is determined by classical 
trajectory simulations for all relevant rotational states 
using Monte Carlo sampling of the initial phase distri- 
butions [2^. Including the resulting populations in the 
laser interaction zone, the final 2D-screen probability dis- 
tribution reads 



P^{y,z)=Y,W,P;{x,y), 



(5) 



0.4 




Fig. 4: Experimental recoil velocity distribution of CN+ ions 
rescaled to have a maximum for a = 1 (points), fitted func- 
tions using two Gaussians (solid line) and two Gaussians with- 
out background (dashed line). 



where the sum runs over all populated states and the 
index i indicates the polarization of the probe. 



D. The experimental observables 



The alignment of a molecule is quantified by the 



expectation value (cos 



with < 



(COS' 



< 



1, where larger values correspond to stronger align- 
ment. While ((cos^0)) = J^-y^-yi^^^"^ ^)'y is not ex- 
perimentally determined, the alignment is character- 
ized through ((cos^02n)) = ^jW^lcos^ 92d)-i, where 
= arctan(z/j/) is the angle between the z-axis of 
the screen plane and the projection of the ion recoil ve- 
locity vector onto the detector plane. Let us remark that 
(cos^ 62 0) -y is compu ted from the final 2D probability 
density equation (5)[ whereas (cos^ 9)^ from the 3D wave 
functions. 

When the linear polarization of the YAG laser is not 
perpendicular to the static electric field (/3 7^ 90°), the 
up/down symmetry of the 2D- images is lost, and an 
asymmetric distribution appears showing a certain de- 
gree of orientation. This up/down asymmetry is experi- 
mentally quantified by the ratio iVup/TVtot, with A^up be- 
ing the amount of ions in the upper part of the screen 
plane, and A'tot the total number of detected ions. In our 
description, they are given by 



and 



A^„, 



Priy, z)dydz, 



(6) 



z>0 



00 /'OO 



PT{y,z)dydz. 



Note that due to the normalization we have that Ntot = 
1. 

Finally, to compare our numerical data with the ex- 
perimental results, we take into account the finite spatial 
width of the YAG and probe pulses. Both laser beams 
have a Gaussian profile with widths of LOy — 36 /im and 
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LUp = 21 /zm, respectively, and are overlapped in time 
and space. The YAG pulse ensures the alignment or 
orientation of the molecules, whereas the probe pulse is 
needed for detection purposes. The degree of alignment 
and the dissociation probability varies with the position 
of a molecule in the interaction volume. Thus, we in- 
tegrate over the overlap region of the probe laser with 
the molecular beam considering the YAG-laser intensity. 
Hereby, we have assumed a linear behavior of the dis- 
sociation efficiency on the probe laser intensity. This is 
an approximation, as recent works about non resonant 
strong field dissociation for other systems have proved a 
/■^-dependence [i^. However, our calculations indicate 
that the orientation and alignment results are not very 
sensitive to a variation of this dependency. We have also 
assumed that the spatial profile of the molecular beam 
is much broader than that of the laser beams, hence, the 
variations of the density of molecules can be neglected. 



IV. COMPUTATIONAL RESULTS 

In this section, we apply the above described approach 
using benzonitrile as a prototype example. Recent exper- 
imental results for this molecule [3l| allow us to directly 
compare them to our numerical studies. The moments 
of inertia are given by Bxm — 1214 MHz, By^ = 1547 
MHz, and Bz^ — 5655 MHz, the electric dipole moment 
is /U = 4.515 D, and the principal moments of polariz- 
ability are axMjf,j = 7.49 A^, ay^iY^i = 13.01 A^, and 
az^z^, =18.64 A3 [mil . 

A. Alignment results 

The experimental fields geometry consists of a weak 
static electric field perpendicular to the screen with 
strength Eg = 286 V cm~^ and a strong alignment laser 
linearly polarized along the Z^-axis (z-axis). For the 
probe pulse we consider the three possible polarizations: 
i) linearly polarized perpendicular to the screen; ii) lin- 
early polarized parallel to the screen; and iii) circularly 
polarized in the x^-plane perpendicular to the screen. 
The computational results shown below include the re- 
coil velocity distribution and the volume effect. 

The theoretical results for the density distribution 
of the molecular beam on the screen are presented 
in [Figure 5| (a), (b) and (c) for Eg = 286 V cm^^ 
and for YAG pulse intensities / = 5 x 10^° W cm^^, 
10^^ W cm~^, and 7 x 10^^ W cm~^, respectively A 
strong confinement of the probability distribution is ob- 
served, the molecules are aligned with their Zjv/-axis 
pointing along the polarization axis of the YAG laser. 
This ID alignment becomes more pronounced as / is in- 
creased. The two ionization channels that characterize 
the velocity distribution for BN, cf. [Figure 4] appear as 
double maxima on the upper and lower hump of the 2D 
images. They become more prominent as / is increased. 




O I 8 I 8 



Fig. 5: The 2D projection of the density distribution of the 
benzonitrile molecular beam for (a) 7 = 5 X 10^° W cm-^ (b) 
10" W cm-2 and (c) 7 x lO" W cybT'^ , Es = 286 V cm"^ 
and /3 = 90°. 
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Fig. 6: Alignment cosines for the molecular beam of be- 
zonitrile as a function of the YAG pulse intensity, Es = 286 
V cm~^ and /3 — 90°: {{cos^ 6)) (circles and dashed line), 
((cos'^ 62d)Y (squared and solid line), {{cos^ 62d)Y (asterisks 
and dotted line), and ((cos^ 620))^ (triangles and dash-dotted 
line) . 



i.e., for a YAG pulse with / = 7 x 10" W cm'^^ while 
they are overlapping for lower intensities. 

We quantify this alignment by the expectation values 
{{cosH)), ((cos2 02Z5))M(cos^2D»^ and {{cosH2D)r, 
which arc presented in [Figurc~6| as a function of the YAG 
intensity / for the whole molecular beam. The indexes 
I, p and c indicate the probe pulse polarization: linearly- 
polarized perpendicular and linearly-polarized parallel 
to the static field and circularly-polarized, respectively. 
Without the presence of any aligning or orienting fields 
the screen image is already biased by the "geometric 
alignment" of the probe pulse. An analytical calcu- 
lation shows for a random isotropic distribution that 
{{cosH)) = 1/3, ((cos2 620))' = 0.5, ((cos2 02d))p = 3/4, 
and ((cos^02d))^ = 5/8. This mild confiment obtained 
when the probe pulse is circularly polarized or linearly 
polarized parallel to the screen is due to the enhanced 
ionization probability for molecules having the C2 axis 
(i.e., the C-CN axis) parallel to the laser polarization. 
If the molecules are exposed only to the weak electric 
field, Es = 286 V cm~^, these field-free values are only 
slightly perturbed ((cos^ 6)) « 0.34, ((cos^ 6*2^))' = 0.51, 
((cos2 6I2d))p = 0.7496, and ((cos^ 620))" = 0.6319 
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Fig. 7: The 2D projection of the density distribution of the 
benzonitrile molecular beam for (a) / = 10", and (b) 7x lO" 
W cm-2, Es = 286 V cm~^ and /3 = 135°. 



These four expectation values qualitatively all show 
a similar evolution as the YAG pulse intensity is var- 
ied: a steep increase followed by a plateau-like behav- 
ior which indicates a saturation of the degree of align- 
ment. This dependence of the alignment on the YAG 
pulse intensity nicely reproduces the experimental behav- 
ior observed for a molecular beam of iodobenzene (IB) 
[29| . Note that the differences due to the polariza- 
tions of the probe pulse are noticeable for low intensities, 
whereas for strong alignment fields they all converge to 
the same asymptotic limit. For I = 7 x 10^^ W cm~^, 
we obtain {{cos^ 9)) = 0.946, {{005^620))^ = 0.972, 
((cos2 6'2d))p = 0.973, and {{cos^ 620}}" = 0.974. The 
value obtained experimentally, 0.89, for ((cos6'2d))p [31| 
is somewhat lower which we ascribe to contaminant ions, 
like C2HJ (with the same mass-to-charge ratio as CN+) 
in the images. Ions like C2ll^ do not have a strong angu- 
lar confinement and therefore reduced the apparent de- 
gree of alignment in the image. For IB, a molecule which 
is expected to attain an alignment degree similar to that 
of BN because of a similar polarizabilty tensor, the recoil 
ion used, I"*", is more clean since there are no contami- 
nant ions at the mass of 127. As a result, the observed 
degree of ahgnment is as high as 0.97 for ((cos02_d))' HI] 
in good agreement with the theoretical predictions. 

For completeness, we have also considered a thermal 
sample of BN at 1 K and I = 7 x 10" W cm'^, ob- 
taining ((cos^ 02D)y ~ 0.961, which agrees well with the 
experimental value for IB [48| . 



B. Orientation results 

The polarization axis of the YAG-laser is now rotated 
about the l^-axis, and forms an angle /3 {/3 ^ 90°) with 
the weak electrostatic field perpendicular to the screen, 



286 V cm- 



The orientation ratio Nup/Ntot 



equation (6) is derived for a circularly polarized probe 



The theoretical results for the projected density dis- 
tribution of the molecular beam are presented in Figs. 
[Tlja) and (b) for two intensities and /3 = 135°. As a con- 
sequence of the rotation of the YAG polarization axis, 
the up/down symmetries of the wave functions is lost. 
The ID orientation becomes more pronounced as / is in- 
creased. For /3 — 135°, the projection of CN+ ions on the 
2D detector overlaps the two Coulomb explosion channels 
to an extent that they cannot be discerned. 

The theoretical results for the orientation Nup/Ntot are 
presented in |Figurc~8| as a function of the tilt angle /3 and 
for three intensities. The main feature of the Nup/Ntot is 
that it is almost independent of the alignment laser in- 
tensity and of (3. For these three intensities, all the popu- 
lated states are within the pendular regime. Thus, if the 
volume effect was neglected, the ratio Nup/Ntot would be 
independent of / for this pendular regime of laser inten- 
sities. By taking into account the spatial distribution of 
the probe and YAG beams, the molecules show the small- 
est orientation for / = 10^^ W cm~^. This is explained 
by the larger contribution of lower intensities, for which 
the molecules are not yet in the pendular regime, to the 
volume effect integral for / = 10^^ W cm^^, and as a 
consequence Nup/Ntot is reduced. For / = 5 x 10^"'^ and 
7 X 10^^ W cm~^, the contribution of lower intensities 
is not strong enough to cause any noticeable difference. 
Note that for IB, it was also experimentally found that 
Nup/Ntot is independent of the YAG pulse intensity [29| . 
For p = 90°, the molecular beam does not show any 



0.64 
0.62 
0.6 
0.58 
0.56 




0.54 



0.52 



pulse, including the velocity distribution and the volume 
effect and using the diabatic model described above. 



90 95 100 105 110 115 120 125 130 135 

Fig. 8: The theoretical orientation ratio Nup/Ntot as a func- 
tion of the angle /3 for the molecular beam of benzonitrile, 
YAG pulses with / = 10^^ (squares and solid line), 5 x 10^^ 
(circles and dashed line) and 7 x 10^^ W cm~^ (triangles and 
dash-dotted hue) and a static field with Es = 286 V cm~^. 

orientation and it holds Nup/Ntot ~ 0.5. For a certain 
YAG pulse intensity, Nup/Ntot quickly increases as /3 is 
increased from 90 and reaches a plateau-like behavior for 
the rest of the field configurations. For /3 = 135° and / = 
7 X 10" W cm-2, we obtain Nup/Ntot = 0.636 which is 
smaller than the experimental result of 0.71 [3l|. Thus, 
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our theoretical model underestimates the orientation ob- 
tained in the experiment. 

An important discrepancy between theoretical and ex- 
perimental results is the variation of Ny^Ntot for /3 close 
to 90°. For a molecular beam of IB [2^, it was exper- 
imentally shown that Nup/Ntot smoothly increases (de- 
creases) as /3 is increased (decreased) from 90°, reaching 
for P < 60° and P > 120° the plateau. In contrast, 
the theoretical Nup/Ntot sharply increases, and its (3- 
independent value is already achieved for /? > 91°. In the 
pendular regime, the strong laser interaction pairs states 
into quasidegenerate doublets. The electrostatic field can 
induce a strong coupling between these levels giving rise 
to a large orientation if the energy gap is small enough. 
For certain intensities, some excited levels may not show 
such a narrow energy gap to obtain a significant orienta- 
tion, but, their relative weight within the molecular beam 
is so small that their contribution to the final result is not 
relevant. For a linear molecule in combined fields j49l |. 
such a sharp rise was predicted for the expectation value 
(cosfl) as the static field strength is increased, i.e., in- 
creasing or decreasing f] from 90° in our case. On Figure 
2 of Friedrich and Herschbach Ref. [22| this effect is ob- 
tained for the interaction with the laser field being 25000 
times larger than the one with the static field. We could 
perform a similar comparison and for BN, I = 7 x 10^^ 
W cm-2, Es = 286 V cm"! and /3 = 95°, the interaction 
of the YAG pulse, 2TTla^^ /c, is 43600 times larger than 
the coupling with the static field, fiEs cos p. 

Let us remark, that our diabatic model for the popula- 
tion transfer is equivalent to considering a field configu- 
ration with a linearly polarized YAG pulse being parallel 
to an electrostatic field with strength Eg cos /3. This field 
geometry is neglecting the component of the electric field 
Es sin /? perpendicular to the YAG polarization, which 
is responsible for breaking the azimuthal symmetry and 
causing M to stop being a good quantum number. This 
approximation can be done because the interaction due 
to the static field is sufficiently weak. We have performed 
the full calculation considering the simplified case of par- 
allel fields, the static one with strength Eg = 286 cos 135° 
V cm~^ and / = 7 x 10^^ W cm~^, obtaining for the ori- 
entation the value Nup/Ntot = 0.636. 

To illustrate the necessity of the diabatic model, we 
build up the molecular ensemble by successively adding 
states according to their weight. Thus, our initial en- 
semble contains only the ground state, and for the sec- 
ond one, we add the second and third most populated 
levels. The relative weights of the states within these 
molecular beams are the same as in the full molecular 
ensemble. [Figure 9] presents the orientation Nup/Ntot as 
a function of the percentage of states included in the en- 
semble for / = 7 X 10" W cm-2 and /3 = 135°. The ratio 
Nup/Ntot has been computed with the diabatic model for 
the population transfer and a fully adiabatic picture, i.e., 
all avoided crossings arc assumed to be passed adiabat- 
ically. The first three populated states, OqoO, IqiI and 
loi — 1, are perfectly oriented: N^p/Ntot = 1- For the sets 
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Population percentage 

Fig. 9: The theoretical orientation ratio Nup/Ntot as a 
function of the population of the molecular beam of ben- 
zonitrile, for the field configuration J = 7 x 10^^ W cm~^, 
Es = 286 V cm~^, and /3 — 135°, computed with diabatic 
model for the population transfer (squares and solid line) and 
a fully adiabatic picture (circles and dashed line). 

including the levels adding 20% of the total population, 
the diabatic and adiabatic results agree because these 
states do not suffer any avoided crossings. As highly 
excited levels are added to the molecular beam, the dif- 
ference between both calculations becomes more evident 
due to the presence of the diabatic avoided crossings. 
They both show a zig-zag decreasing trend as the popu- 
lation is increased. The ratio Nup/Ntot computed under 
the adiabaticity assumption is smaller than the diabatic 
result as more states arc added. Taking into account all 
states populated in the molecular beam, the adiabatic 
model does not yield any appreciable amount of orienta- 
tion: we obtain Nup/Ntot = 0.53 which strongly underes- 
timates the experimental result. In addition wc have ob- 
served that the adiabatic orientation ratio Nup/Ntot is no 
longer /^-independent. This confirms that non-adiabatic 
crossings play a crucial role for understanding the strong 
orientation observed in the experiment. 

Finally, it is worth noting for the field configuration 
/ = 7 X 10" W cm"2 and (3 = 135°, a thermal sample 
of BN at 1 K shows a weak degree of orientation with 
Nup/Ntot = 0.563. 



V. SUMMARY AND CONCLUSIONS 

We have presented a theoretical model to investigate 
the degree of alignment and orientation of a beam of 
asymmetric top molecules exposed to combined electro- 
static and non-resonant linearly polarized laser fields. 
Our analysis combines the field-dressed wave functions 
with the experimental distribution of the populated 
quantum states. As a first step, we solve the three- 
dimensional Schrodingcr equation within the rigid ro- 
tor approximation. For a certain field geometry, we 
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treat each irreducible representation independently, by 
expanding the wave function in a basis that respects the 
corresponding symmetries. Since the dc electric field 
strength is very weak, we consider M as being almost 
conserved, and a diabatic model is introduced to approx- 
imate the population transfer trough the avoided cross- 
ings as the YAG pulse intensity is increased. The 2D 
projection of a wave function is derived by using the de- 
tection selectivity factors of the probe pulse, the velocity 
distribution of the detected ions, and a volume effect av- 
erage. 

This theoretical model has been checked by comparing 
the numerical and experimental results for benzonitrile. 
The molecular beam has been simulated using the popu- 
lation of each quantum state numerically obtained from 
deflection profiles [2^. For several field configurations, 
we have performed a detailed study of the directional 
properties of the molecular mixture. In particular, we 
have explored the degrees of alignment and orientation 
as the YAG pulse intensity and the angle between both 
fields are varied. For perpendicular fields, a good agree- 
ment between the computational and experimental re- 
sults is obtained. Let us remark that we do not take into 
account the background of unwanted ions which contam- 
inate the velocity mapping images for BN and reduce 
the experimental degree of alignment [3l|. Hence, a bet- 
ter agreement could be achieved for other systems with 
cleaner Coulomb explosion imaging channels, like IB. 

Regarding the orientation results, we have shown that 
the assumption of a fully adiabatic dynamics is incor- 
rect for the prototypical experiment and predicts a non- 
oriented molecular beam - which is not in agreement with 
the experimental results. Indeed, we have proven that 
the degree of orientation does not adiabatically follow 
the time envelope of the YAG laser with a FWHM= 10 
ns in the experiment. By employing a simple diabatic 
model, the experimental results for orientation could be 
reproduced reasonably well. Based on the comparison 
with the experimental measures, the important impact of 
the diabaticity on the field-dressed molecular dynamics 
is hereby demonstrated. However, notice that our model 
does not produce the smooth /3-dependence of Nup/Ntot 
that was experimentally obtained for IB [29| . 

This theoretical model is based on several approxima- 
tions that may be the source of discrepancies with the 
experimental results. The molecules are exposed to an 
alignment laser pulse with a certain time profile, but 
we are performing a time-independent description of the 
field-dressed rotational dynamics. We compensate this 
deficiency by using a simple partially diabatic approach 



to model the populations transfer. However, a general 
statement on the character of the avoided crossings can 
not be made, and some of them could not be classified 
clearly by our adiabaticity criterion. The semiclassical 
calculations of the molecular trajectories through the de- 
fiector have been derived using a fully adiabatic picture, 
whereas it has been theoretically shown that such an ap- 
proximation may be incorrect [39j. These assumptions 
for both processes might not be fully satisfied. A full 
time-dependent description of the process would properly 
treat the avoided crossing, but it is computationally very 
challenging. Possible distortion effects due to the strong 
laser pulse are also not taken into account. Furthermore, 
we are working within the axial recoil approximation ne- 
glecting the interactions between the ionic fragments on 
their way to the detector. Analogously, collisions and 
interactions between the molecules within the molecular 
beam [50| were not considered. The geometric alignment 
due to the strong-field ionization has been assumed to be 
a single-photon absorption, but it is a multiphoton pro- 
cess [46|. Nuclear hyperfine structure has been neglected. 
Whereas, it can be of the same order of magnitude as the 
interaction with the weak static field (i.e., for IB [sij), 
this is clearly not the case for BN 

Certainly, it would be interesting to perform a com- 
parison of these computational results with experimen- 
tal data obtained for other molecular species either with 
a smaller number of populated states or with a cleaner 
imaging signal. A rather natural extension of the present 
work would be to go beyond any of the approxima- 
tions described above. A more sensitive criterion for the 
avoided crossing treatment, a time-dependent descrip- 
tion, or the inclusion of the hyperfine interaction should 
improve these theoretical results. 
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